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We investigate the stability of few-electron quantum phases in vertically coupled quantum dots 
under a magnetic field of arbitrary strength and direction. The orbital and spin stability diagrams 
of realistic devices containing up to five electrons, from strong to weak inter-dot coupling, is de- 
termined. Correlation effects and realistic sample geometries are fully taken into account within 
the Full Configuration Interaction method. In general, the magnetic field drives the system into a 
strongly correlated regime by modulating the single-particle gaps. In coupled quantum dots difi'erent 
components of the field, either parallel or perpendicular to the tunneling direction, affect single-dot 
orbitals and tunneling energy, respectively. Therefore, the stability of the quantum phases is re- 
lated to different correlation mechanisms, depending on the field direction. Comparison of exact 
diagonalization results with simple models allows to identify the specific role of correlations. 
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I. INTRODUCTION 

Due to three-dimensional confinement on length scales 
comparable to the De Broglie wavelength, the electronic 
properties of semiconductor quantum dots (QDs) show 
several similarities with those of atoms, the most sig- 
nificant ones being their discrete energy spectrum and 
the resulting shell structurei^ Fine structure due to 
the exchange interaction (Hund's rule)^*^*^ and Kondo 
physics^ have also been predicted and demonstrated. 
QDs are therefore regarded as artificial atoms^'^'^'iS^ 
These systems have stimulated the investigation of the 
fundamentals of few body physics in semiconductors, 
since the number of carriers, electrons or holes, in QDs 
can be controlled very accurately, and almost all relevant 
parameters influencing their strongly correlated states, 
such as confinement potential or the coupling with static 
magnetic fields, can be tailored in the experiments, driv- 
ing the system between very different regimes. 

The analogy between artificial and natural atoms is ex- 
tended to artificial molecules (AMs) by realizing tunnel- 
coupled devices jiSiiii^ii^iiSiii The inter-dot tunnelling 
introduces a new energy scale which may be comparable 
to other energy scales, namely the single-particle (SP) 
confinement energy, the carrier-carrier interaction, and 
the magnetic confinement energy. In AMs it is possible 
to tunc the coupling among QDs basically at will, thus 
exploring different molecular bonding regimespS- this is 
an intriguing option with respect to natural molecules, 
since for the latter the inter-nuclear coupling is almost 
fixed by the balance between nuclear repulsion and the 
electrostatic attraction mediated by valence electrons. 

There is a number of different techniques to fabricate 
AMs. Laterally coupled AMo^^i^°i^^i^^i^^ are obtained 
by creating an electrostatic confinement in a semicon- 
ductor heterostructurc, such as a doped heterojunction, 
by means of the photo-lithographic patterning of metal- 
lic gates, deposited on the surface of the heterostruc- 
turc. Another way to realize AMs is by the Stransky- 



Krastanov mechanism leading to the formation of self- 
assembled QDs with nm-scale confinement with similar 
sizes and regular shapes. In a stack of several layers, the 
formation of QDs in the top layer, coupled through tun- 
neling with the ones in the underlying layer, has been 
recently demonstratediS^iS^ 

In the following we shall investigate electronic prop- 
erties of vertically coupled AMs realized starting from 
a triple barrier heterostructurc to form a mesa pil- 
lar, possibly using a combination of dry and wet 
etching, 141^8,26,27,28,29,30,31 Typically, the effective diam- 
eter for electron conduction is smaller than the geometri- 
cal diameter of the top contact. For example, with a top 
contact « 0.5 /im diameter, it is estimated that the dot 
diameter is ~ 0.1 fim when the dot contains a few elec- 
trons. Still, in this class of devices the depletion potential 
in the lateral direction realizes a much softer potential 
than in the growth (vertical) direction, where electro- 
static confinement is induced by the band mismatch of 
the layers in the heterostructurc. In this vertical geome- 
try transport occurs perpendicularly to the plane of the 
dots, in response to an applied voltage between source 
and drain contacts. The strength of the coupling can 
be modulated by the width of the barrier separating the 
QDs. In addition, a circular Schottky gate is typically 
placed around the body of the pillar to control the charg- 
ing of the QDs. 

In AMs, as well as in single QDs, electronic states can 
be manipulated by a vertical magnetic field B± , perpen- 
dicular to the plane of the QDs (parallel to the growth 
direction), which drives the system from a low-correlation 
(low-field) regime to a strongly correlated (high-field) 
regime by varying the SP splittings.^ Recently, non triv- 
ial transitions between different quantum states induced 
by B± and/or by varying the inter-dot barrier width if, 
have been predicted and demonstrated in AMs with dif- 
ferent coupling regimes liSiiSi^SiS^ 

While this vertical field configuration has 
been widely studied both theoretically and 
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few works have been devoted to the in-plane field 
configuration(21iiiiiSiSiiS which is more complicated 
to simulate and less intuitive, principally because of 
the lack of the cylindrical symmetry and of analytical 
solutions for the SP states. Furthermore, since in this 
configuration angular momentum is not a good quantum 
number, the computational effort increases very rapidly 
with the electron number (see Sec. II B). 

On the other hand, in vertically coupled QDs an in- 
plane magnetic field _B|| modulates the tunneling en- 
ergy, and one may hope to induce transitions between 
few-particle states, reminiscent of those that one can 
obtain by varying the tunneling energy via the barrier 
potential^ In this way, different coupling regimes would 
be investigated within just one sample, which is of course 
an advantage, since it is almost impossible to growth 
samples completely identical except for the barrier width. 
It should be stressed that when the tunneling is modu- 
lated by an in-plane magnetic field in vertically coupled 
QDs, the other energy scales, in particular the Coulomb 
interaction, are only weakly affected by the field. Due 
to its importance to the QD-based implementation of 
quantum-information processingji^ such field configura- 
tion has been already investigated in the two-electron 
case^ whereas less attention has been devoted to the 
case where more than two carriers are present, in par- 
ticular when B has both an in-plane and a perpendicu- 
lar component, which is relevant for the transport spec- 
troscopy of such samples. In this case we expect an in- 
teresting physics: in fact both in-plane (due to B±) and 
vertical (due to _B||) correlations are involved in the tran- 
sitions between quantum states. 

In this paper we theoretically investigate the stabil- 
ity of field-induced quantum states of few electrons in 
AMs. The main focus is on correlation effects, largely 
enhanced by the field, which determine different mecha- 
nisms (which may interfere between each other) driving 
the transition between different ground states, depend- 
ing on the field direction and inter-dot coupling. In our 
analysis we consider three realistic samples with differ- 
ent inter-dot barrier, with the aim of investigating the 
phenomenology occurring between the two opposite lim- 
its of strong coupling and molecular dissociation. Our 
numerical approach is based on a real-space description 
of SP states, which takes into account the complexity of 
real samples. Because of the competition between kinetic 
energy, Coulomb interaction, and Zeeman energy, it is of- 
ten necessary to treat exactly (in contrast, e.g, to mean- 
field methods) the few-body Hamiltonian of this kind of 
samples. Our method of choice to include carrier-carrier 
Coulomb interaction is the Full Configuration Interac- 
tion (FCI) approach, which proved to be accurate and 
reliablejiiiSiSiS although limited to few electrons. In 
order to stress correlation effects, we contrast our nu- 
merical results with SP and Hartree-Fock predictions of 
the stable quantum phases. We will see that correlation 
effects play a crucial role in determining the electronic 



^2i^#?^rtId*'c*?-??-^t*^tem, and that schematic pictures ne- 
glecting these phenomena do not lead to correct predic- 
tions. 

The paper is organized as follows: In Sec. II we illus- 
trate our numerical approach to the calculation of SP and 
interacting states; in Sec. Ill we discuss the results con- 
cerning different samples in a magnetic field of arbitrary 
direction; in Sec. IV we summarize our findings. 



II. THEORETICAL MODEL 

We consider N interacting electrons in an AM struc- 
ture. Since we are interested in weakly confined devices, 
the energy region of our concern is relatively close to 
the semiconductor band gap and we may use the single- 
band approximation for the conduction band. In this 
framework the carriers are described by the effective- 
mass Hamiltonian 



N 



^(-TO + MA(rO)%nr,:) 



■E 



9 MB 



K r,; - To 



B-S. (I) 



Here to*, k, and g* are the effective mass, dielectric con- 
stant, and 5-factor, respectively, is the Bohr magne- 
ton, S is the total spin, A(r) is the vector potential at 
position r, and B is the magnetic field, with B = V x A. 
Equation Q neglects non-parabolicity effects, but it in- 
cludes otherwise in the confinement potential V{r) the 
full three-dimensional nature of quantum states in real- 
istic samples, such as layer width, tunneling, and finite 
band offsets. 

In the following the potential V^(r) describes two iden- 
tical vertically coupled disk-shaped QDs. Since in this 
kind of samples the confinement is much tighter in the 
growth direction z than in the xy plane, we separate the 
potential as V{r) — V{x,y) + V{z), where V{z) repre- 
sents two identical square quantum wells of width Lw, 
separated by a barrier of width if, and conduction band 
mismatch Vq . We perform the usual choice of a parabolic 
in-plane confinement of natural frequency loq: 



V{x,y) = ^m*cul{x^ +y^). 



(2) 



Adequacy of Eq. Q has been demonstrated by both the- 
oretical calculations^^"^'' and far infra-red spectroscopy 
experiments.^^'^- Note, however, that our numerical ap- 
proach does not assume any symmetry and may treat 
different confinement potentials; in particular, the vector 
potential A(r) is not limited to describe z-directed field. 
We suppose that the active region of the AM is suffi- 
ciently well separated from the top and substrate con- 
tacts that we can reasonably neglect interactions with 
the leads in our modelfiS, The energy splitting Tiuiq be- 
tween SP eigenvalues of the in-plane confinement,"'^'^ given 
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by Eq. is an input parameter of the model. Such 
quantity mimics the effects of electrode screening and 
interaction with the environment, and in real charging 
experiments it is modified by the Schottky gate.^^ Here, 
for the sake of simplicity, Hloq is kept fixed as N is varied. 
A qualitative discussion of this issue is postponed to the 
end of Sec. |llll 

Schematically, our algorithm is the following. First, we 
numerically calculate the electron SP states, mapping the 
SP Hamiltonian on a real-space grid. Then, the SP or- 
bitals thus obtained are used to evaluate Coulomb matrix 
elements and finally to represent the interacting Hamilto- 
nian on the basis of Slater determinants (SDs), according 
to the FCI approach. These steps are described below in 
detail. 



A. Single-particle states 

The SP energies and wave functions 0c[(r) are ob- 
tained from the numerical solution of the eigenvalue prob- 
lem associated to the Hamiltonian 
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which is the SP term appearing in the first line of Eq. Q • 
We include in Eq. |j2Jl a magnetic field of arbitrary direc- 
tion via the vector potential A = i (B x r) (symmetric 
gauge). Specifically, we find the eigenstates of Eq. Q 
by mapping it on a real-space grid of A'g^id = Y[i=i 
points, identified by the grid vectors r.^ — X^fc^il'^i^ ^ 
7Vfc/2)Afcefc, with A*^ = l,...,Nk and §1,2,3 = x,y,z. 
The resulting finite-difference equation can be rewritten 
in terms of a discrete eigenvalue problem, as described in 
detail elsewhereiS This results in the diagonalization of 
a large, sparse matrix which is performed by the Lanczos 
methodi^ 



B. Few- particle states 

Once SP states are obtained, we proceed including the 
Coulomb interactions between charge carriers. The full 
many-body Hamiltonian H can be rephrased, by means 
of the second-quantization formalism, as 

H = Ho + Hc + Hz, (4) 

where Hq is the single particle Hamiltonian 
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Hn is the Coulomb term 
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and Hz is the Zeeman Hamiltonian 
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Here, c[^^ {caa) creates (destroys) an electron in the spin- 
orbital 0Q(T(r, s) — (l)a{v)xa{s), where Xct(s) is the spinor 
wave function, and S is the total spin vector. In the 
above equations, A^sp is the number of SP states that 
are taken into account. 

The Coulomb matrix elements, given by 



drdr', (8) 



are calculated by numerically integrating the following 
expression 
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where 3'a/3(r) = (/>* (r)0/3(r), and $a/3(k) = J^[$a/3(r)] is 
its Fourier transform. 

Since typical SP energy spacings computed from the 
Schrodinger equation associated with Eq. (j2Jl are in the 
meV range, while characteristic values of the Zeeman 
term Hz are typically two orders of magnitude smaller 
(e.g., [i^B = 5.79 X 10"^ meV at 1 T), the effect of Hz is 
first neglected in the calculation of the electronic ground 
states (GSs). Then, the effect of Hz, which just lifts the 
(25" -I- l)-fold degeneracy of the total spin 5*, is included 
as a perturbation to the first order in the field. 

Then, our algorithm proceeds forming a basis of SDs, 
filling iVsp spin-orbitals (j)aa{T^, s), calculated numerically 
as described in Sec. II A, with N electrons in all possible 
ways. We assume that the Fock space generated by the 
SDs is approximately complete, namely a generic eigen- 
state of the many-body Hamiltonian H can be expressed 
as a linear combination of SDs. One may reduce the 
computational effort by exploiting the symmetries of the 
system. In fact, by using a suitable combination of the 
SDs which diagonalize symmetry- related operators, the 
Hamiltonian can be written in a block-diagonal form.^ 
The larger the symmetry of the system, the smaller the 
dimensions of the blocks and the time necessary to com- 
plete the full diagonalization. In the present case, an 
in-plane component of the magnetic field removes the 
cylindrical symmetry Cooh, which is often exploited in 
the vertical field arrangement of cylindrically symmetric 
QDs. However, the Hamiltonian still commutes with the 
square total spin S'^ and z-component Sz operators, re- 
spectively (neglecting Hz)- Therefore, the subspaces are 
appropriately labeled by the values of S and the mini- 
mum positive value of Sz consistent with S. Once the 
effective subspace is selected, our code performs a uni- 
tary transformation in order to rewrite H in a block- 
diagonal form. Finally, matrices obtained in this way are 
handled by a Lanczos routine included in the package 
donrodrigc^^ run on a SP3 IBM system. 
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FIG. 1: Fock-Darwin energy levels of samples Bl, B2, B3 
(see text) vs. vertical magnetic field, _Bj_. Solid and dashed 
lines represent energies e„m arising from the lowest S and 
AS levels of the double quantum well confinement along the 
growth direction, respectively. 



III. RESULTS 

In this section we present the calculated field- 
dependent GSs of AMs, covering the regimes from strong 
to weak coupling. The appUcation of the FCI method en- 
sures that our results have comparable accuracy in differ- 
ent regimes. To this end, we have selected three realistic 
samples based on a AlGaAs/GaAs heterostructure, la- 
beled Bl, B2, B3, with inter-dot barriers Lfc = 2.5, 3.2, 
and 4.0 nm, respectively. The sample with intermediate 
barrier shows a full molecular character, while Bl and 
B3 represent the siglc QD and the molecular dissocia- 
tion limits, respectively. 

Except for the barrier width, parameters are common 
to all samples. The lateral confinement is htoo — 3.5 meV, 
the band-offset between barrier and well materials of the 
heterostructure is Vq = 300 meV, each well width of the 
double quantum well potential is Lw = 12 nm. Material 
parameters appropriate to GaAs are m* /rrie — 0.067, 
K = 12.9, g* = —0.44. All parameters refer to a set of 



FIG. 2: Single-particle states of samples Bl, B2, B3 vs. in- 
plane field -B|| at B± = 0. 
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FIG. 3; (Color online) Contour plots of \(f){x,y)f, at a fixed 
value of z inside one well, of selected SP states (see labels 
on top) of sample B3 with an in-plane field -By = (top 
row), -By = 4 T (center row) and -By = 8 T (bottom row). 

lies along the x axis. Symmetry labels refer to the lowest- 
energy SP orbitals of a given representation (cf. bottom panel 
of Fig.Hll. 



5 



sample3i2i2i described in literature and processed at the 
University of Tokyo. 



A. Single-particle states 

SP states in an arbitrary field have been obtained ac- 
cording to Sec. Hm with Ni = N2 = 80, N3 = 128; the 
resulting 819200 x 819200 sparse Hamiltonian matrix is 
diagonalized by means of the Lanczos algorithm. The 
grid parameters have been chosen in order to reproduce 
the lowest analytical zero field (Fock-Darwin) energies 
£nm with an accuracy better then 0.5 meV. 

Before presenting our results for specific samples, let 
us summarize the properties of SP states in a single QD 
with parabolic in-plane confinement and vertical mag- 
netic field B±. In this case SP states are analytical and 
given by the so-called Fock-Darwin states (FD) , with en- 
ergies Enm = fifl{2n+\m\ -I- 1) — (hwc/'2.)m, n = 0, 1, 2, . . . 
and TO = 0, ±1, ±2, . . . being the principal and azimuthal 
quantum numbers, respectively. The effective oscilla- 
tor frequency Vl is defined by = + '-^c/^j and 
LOc = eB/m*c is the cyclotron frequency. A perpendic- 
ular magnetic field splits the degeneracies of the states 
and reduces the energy gap between energy levels with 
increasing angular momentum (see solid lines in Fig.^; 
this favors transitions in the few-electron states in order 
to increase in-plane correlations, as explained in the next 
section. In symmetric AMs, FD levels are replicated, 
rigidly shifted by AsaSj the energy gap between the low- 
est symmetric (S) and anti-symmetric (AS) states arising 
from the double- well potential along the growth direction 
(highest confined states along the growth direction need 
not be considered if Lw is sufficiently small) Figure Q] 
shows the S (solid lines) and AS (dashed lines) replicas 
of FD states for the three samples as a function of B±. 
Note that the S/AS labeling is valid only as far as the 
field is entirely along the z axis, namely the motions in 
the xy plane and along z are completely uncoupled, and 
the splitting Asas does not depend on B± . The symme- 
try group, in this highly symmetric case, is Cooh- 

Figure 12 shows calculated SP levels for Bl, B2, B3 
as a function of the in-plane field . Contrary to the 
vertical-field configuration, the parallel field , which 
we take to lie along the x axis, couples the motion along 
the in-plane y and the vertical z directions, and it re- 
duces the symmetry group of the system from Cooh to 
C2/1. Therefore, SP wave functions lose their well de- 
fined component of the orbital angular momentum as 
well as the S/AS character along z. Besides, the states 
labeled S and AS when = 0, respectively, get closer as 
increases, and tunneling is progressively suppressed. 
Clearly, the tunneling suppression occurs at different 
values of the field for different samples, since the in- 
plane field may significantly affect tunneling only when 

Lul =eB||/TO*c~ AsAs/?i— 

Although the application of strongly reduces the 
symmetry group of the system, it still preserves two resid- 
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FIG. 4: Calculated single-particle energy levels of sample B2 
vs. field, B — i3||X -I- B_i_z, forming a 45° angle with respect 
to the growth direction of the heterostructure. 

ual symmetries: i) the reflection x — > —x with respect to 
the yz plane perpendicular to the field direction ii) the 
rotation of 180° around the x axis. The four irreducible 
representations of the residual group C2h, namely Au, 
Ag, Bu, and Bg^ are associated to the appropriate SP 
states in Fig. |21 together with the alternative labeling at 
zero field. It is possible to understand qualitatively the 
behavior of the SP states with the parallel field by 
taking into account the energy gap between the orbitals 
with the same symmetry; the smaller the gap, the larger 
the "repulsion" between levels. For example, we expect 
strong repulsion between the second and the fourth state, 
and a weaker effect for the first and the fifth state because 
of the larger energy gap between them. 

The effect of the in-plane field on SP states is demon- 
strated in Fig. 131 which shows the SP charge densities 
|(/)(x,?/)p at a fixed value of z inside one well for the 
three lowest levels (their quantum numbers are indicated 
in the figure), for three selected values of -By. The field 
increases the confinement, squeezing the electronic wave 
functions. Moreover, FD orbitals with different angular 
momentum and belonging to different energy shells are 
mixed by , leading to charge density modulation and 
nodal surfaces. 

In presence of a magnetic field of arbitrary direction 
(tilted magnetic field) also the C2h symmetry is broken, 
and we observe simultaneously the reduction of both tun- 
neling and intradot energy gaps between FD states as the 
field is increased 1^ An example of the SP states behavior 
in a magnetic field which is tilted at 9 — 45° with respect 
to the growth direction is reported in Fig. 21 

B. Few- particle states 

In the following, we will investigate the stability of 
few-particles (2 < < 5) quantum phases as a func- 
tion of an applied field with both in-plane and vertical 
components. In AMs, due to tunneling and the ensuing 
formation of bonding and anti-bonding levels, the num- 
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FIG. 5: A'' = 2 stability diagram in the {B^^,B±) plane for 
Bl, B2, and B3, as predicted by the FCI calculations. Gray 
squares: 5 = 1 GS; black circles: S — GS. Solid lines 
are guides to the eye. The dashed line in B2 indicates the 
predicted SSD singlet/triplet boundary. Insets: main compo- 
nents of the ground state wave functions (see text) . 



ber of SP orbitals which ensures convergence of the FCI 
calculation depends on the splitting, AgAS, and on the 
number of charge carriers, N , we take into account. Re- 
sults presented in this work are obtained using Ngp — 20 
for 2 < < 4. In order to limit the computational 
effort, which increases very rapidly with the number of 
electrons, we reduce to A^sp = 15 the set of SP states 
for N = 5. With these values, the calculated SP states, 
which at arbitrary values of the field are in general non 
degenerate, turn continuously into close zero field degen- 
erate shells. The convergence of the FCI calculations at 
zero field with respect to the selected SP states has been 
carefully investigated in Ref. |5^ 

In order to stress the role of correlations, we will con- 
trast the exact correlated states, obtained by the FCI 
method, with i) the quantum phases predicted in a SP 
picture, i.e., filling SP states according to the Aufbau 
principle, and ii) a Hartree-Fock-like approach. 

In the SP picture the few-electron non-interacting GS 
energy is 
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(10) 



the occupation number of the a-th SP orbital with spin 
a. The second model is reminiscent of the Hartree-Fock 
approach: In this approximation, we calculate the energy 
of an electronic configuration as the expectation value 
of the interacting Hamiltonian on the most weighted SD 
singled out from the FCI expansion of the correlated wave 
function. We call this single SD (SSD) approximation. 
The SSD energy is given by 
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Here the interaction between charge carriers is described 
by the integrals 



Ual3 — 
Kafj = 
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where is the SP energy of the a-th orbital and is 



Wi - r2| 
X 0a(r2)'/>/3(ri)dridr2, 

where Uap is the direct Coulomb integral, accounting for 
the repulsion between two electrons occupying orbitals a 
and /?, and Kap is the exchange integral, giving the ex- 
change interaction between electrons with parallel spins. 
Uaf3 and Kaf3 correspond to Vappa and Va/3a/3 in Eq. JHJ, 
respectively. Note that in both SP and SSD models the 
GS is given by just one configuration, while, in general, 
the true GS is a linear combination of different electronic 
SDs, due to correlation. We also neglect the Zeeman 
coupling to the field at this level. 

Below we shall discuss the GS of few electrons as a 
function of the magnetic field intensity and direction. 
In several cases, we will indicate schematically (see fol- 
lowing figures) the GS configuration, as predicted either 
by the FCI method or by the SSD calculation, in terms 
of arrows, pointing either upwards or downwards, filling 
in either left or right boxes (with respect to a diagonal 
line). The arrows stand for electron spin, while super- 
posed boxes indicate SP orbitals of increasing energy. In 
the vertical-field configuration, left and right boxes indi- 
cate S and AS levels, respectively (right boxes are not 
shown if no AS states are occupied). If a finite iJy com- 
ponent is present, the boxes indicate the orbitals which 
evolve in a continuous manner from those at -B|| = as 
is switched on (see Fig.|21). For example, the lowest- 
energy S (AS) s and p levels at = evolve continu- 
ously into Ag and Au {B^ and Ag) levels, respectively, 
as increases, while remaining identified pictorially by 
the same boxes. For FCI calculation, we show only the 
most weighted configurations. In some cases FCI and 
SSD methods predict the same stable phase, in which 
case we use dotted boxes. When the two methods dis- 
agree, we use solid boxes and indicate the exact (FCI) 
result only. Typically, for a given sample and number 
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of carriers, the set of stable GSs predicted by the SSD 
method is a subset of the FCI stable phases. 

In all stability diagrams shown below, calculated GSs 
are indicated with dots, and lines are only a guide to 
the eye through the calculated points. Due to numeri- 
cal limitations, mainly related to the lack of cylindrical 
symmetry which requires numerical calculation of the SP 
states, we have limited our calculations to a relatively 
coarse grid of points in the {B±, plane. Although we 
believe that we have determined most of the GSs which 
are stable in the different regimes, it is possible that we 
miss some phases which are stable in a small range of 
fields. 

We first consider the two-electron casejSS, In a SP pic- 
ture, the lowest SP state is doubly filled in a singlet 
{S = 0) configuration which is stable in all field regimes, 
since no level crossing occurs in the lowest SP state (see 
Figs.nand|2Jl. Figure|Slshows, instead, that FCI predicts 
a singlet /triplet transition at finite fields, both along the 
B± and By axes. For B2 the SSD phase boundary is also 
reported in Fig. [S] (no qualitative differences are expected 
for the other samples). 

Let us focus for the moment on the SSD stability di- 
agram of Fig. |31 At low magnetic field the ground state 
has a singlet character, with two electrons sitting on the 
same orbital. As the field is increased, we observe a sin- 
glet/triplet transition, which is different in character de- 
pending whether B± or _B|| is varieduSS In the former case, 
only S levels s and p (equivalent to Ag, levels, respec- 
tively) are involved: the triplet configuration minimizes 
the Coulomb interaction by promoting an electron from 
a s to a p orbital, which is more spatially delocalized, 
i.e., Uss > Usp. In this sense, only in-plane degrees of 
freedom play a role in the _B^-induced transition. The 
"Hartree" energy gain Usp — Uss compensates the cost in 
term of SP energy due to the s p orbital promotion. 
Moreover, the exchange interaction energy gain, given 
by the exchange integral Ksp, favors the "ferromagnetic" 
triplet configuration (S = 1) with respect to the "antifer- 
romagnetic" singlet, when the two electrons sit on the s 
and p orbitals, respectively. Therefore, the electrostatic 
energy of the triplet configuration « Usp — Ksp is further 
reduced with respect to the singlet one, « Usp- 

In the in-plane configuration, instead, the tunneling 
energy separating the lowest-energy Ag and Bu orbitals 
is very effectively reduced by (see Fig. |21l. Both or- 
bitals are therefore involved in the formation of singlet 
and triplet states. In the SSD approximation singlet and 
triplet states, the former made of a doubly occupied Ag 
orbital and the latter of two parallel-spin electrons sitting 
on the Ag and Bu orbitals, respectively, tend to the same 
orbital energy. However, the small exchange interaction 
Ka^b^ and the Zeeman term of the S = 1 state, both 
included in the SSD calculation, weakly favor the triplet 
when the tunneling energy is reduced, above a critical 
value of . 

Although the SSD scheme is able to describe, on a 
qualitative level, the singlet/triplet stability regions, ne- 



glecting correlation effects in the SSD scheme strongly 
underestimates the stability region of the singlet state, 
favoring the polarized phase: in fact, the description of 
the singlet state in terms of two electrons sitting on a 
Ag orbital is poor especially at high field, where an in- 
creasing contribution from the B^ orbital gives rise to a 
spatial correlation along the growth direction, i.e., with 
the two electrons sitting on opposite QDs, as we have 
shown in a previous work^^*Si The vertical correlation 
lowers the energy of the singlet state, making it stable in 
a larger iJy-region, although the Zeeman coupling favors 
the triplet configuration. Indeed, the stability of the sin- 
glet state in FCI may be interpreted as resulting from a 
"superexchange" effect The FCI singlet configuration 
can be better understood after a unitary transformation 
of the orbital basis to states localized on either one or 
the other dot. By schematizing the AM as a two-site 
lattice at half filling, the appropriate conceptual frame- 
work, when the QDs are weakly coupled, is that of the 
Hubbard model)S2i2£ where the system has a tendency 
to "antiferromagnetism" , since two electrons in a singlet 
state can lower their energy by virtually hopping from 
one to the other dot, while the same process is prohib- 
ited for a triplet state by Pauli blocking. When tunneling 
energy does not pay off for antiparallel alignment, the 
triplet phase becomes favorable. 

On the other hand, a SSD scheme would neglect also 
the correlation regarding the triplet state when both B± 
and -B|| are present, namely the central region in the di- 
agrams of Fig. |S1 Here the triplet state is given by the 
sum of mainly two components, the first one involving Ag 
and Au orbitals, and the second one involving Ag and _B„ 
orbitals, with coefficients depending on the field strength 
and direction. In this case the effects induced by both 
_B|| and B± compete in a non-trivial manner, as well as 
vertical and in-plane correlations^ 

Finally, we note that the stability region of the singlet 
state with respect to -B|| depends on the inter-dot cou- 
pling, and it decreases going from Bl to B3. In fact, as 
we have seen before, the in-plane field affects significantly 
the tunneling only when wjj = ei?y/m*c is comparable 
with AsAs/^"^ Consistently with the above interpreta- 
tion of the siglet /triplet transition, the critical value of 
B± does not depend on tunneling strength. 

We next consider the = 3 case. In the SP non- 
interacting scheme (Fig. |SJ) only two phases are present, 
both with S=l/2. The field B^ at which the two GS en- 
ergies cross is given by the intersection of the S s (= Ag) 
andp (= Au) levels (Fig.QJ. Figure [7| shows the stability 
diagram obtained from the FCI calculation; for B2 the 
phase boundary obtained from the SSD approach is also 
shown (dashed line). Although the SSD approximation 
seems to describe reasonably well the dynamics along the 
B± axis, where GS configurations coincide with the FCI 
ones, except a slight shift of the phase boundary at lower 
fields with respect to FCI, this scheme fails completely at 
finite i?y. Indeed, SSD severely underestimates the sta- 
bility region of the unpolarized phase {S=l/2) with re- 
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FIG. 6: = 3 stability diagram in the (B\\ , B±) plane for Bl, 
B2, and B3, obtained within the SP non-interacting scheme. 




If""" ' 





.3=3/2 


B2 




3\0 n\B: 




FIG. 7: A'' = 3 stability diagram in the {B^^,B±) plane for 
Bl, B2, and B3, as predicted by the FCI calculations. Gray 
squares: 5 — 3/2 GS; black circles: S = 1/2 GS. Solid lines 
are guides to the eye. The dashed line for B2 is the SSD 
prediction. Insets: main components of the GS wave functions 
(see text). 



spect to the polarized one (5=3/2). This fact is mainly 
connected to correlation effects, as we shall see in the 
following. 

Let us stick for a moment to the vertical field configu- 
ration. As already for TV = 2, also for N — 3 increasing 
B± determines an enhancement of the in-plane corre- 
lation, leaving unaffected the motion along the growth 
direction. For Bl and B2, electrons in the S orbitals 
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FIG. 8: (a) Lowest energy levels for N = 3 vs. in-plane field 
By at B±—0, for sample B3. (b) Coulomb integrals UAgAg 
and Ub^b^ vs. in-plane field -By at B±=0, for the same sam- 
ple. The Coulomb integrals Uaa are defined in Eq. (I12II . 



undergo a number of transitions where p and d levels are 
successively populated in order to minimize the Coulomb 
repulsion. Eventually, the so called maximum density 
droplet (MDD))72i7?i74,7i? i.e., the densest spin-polarized 
configuration possible, is reached. For B3 this phase is 
not observed at = in the range of B± explored 
(Fig. 13 . 

If -B|| is small or zero, electrons are prevented from 
occupying AS orbitals by the kinetic energy cost Asas- 
As is increased, however, Asas decreases. The range 
of fields where AS occupation takes place depends, of 
course, on the zero field tunneling energy, as it is easily 
seen comparing the GS character at low of the B3 
sample, where electrons occupy both S and AS states, 
with that of Bl and B2 samples, where AS states are 
not occupied. From this point of view, increasing 
is analogous to decreasing Lf,. We will find a similar 
behavior for TV = 4, 5. 

On the other hand, in addition to the reduction of the 
tunneling splitting Asas, a finite value of _B|| determines 
a redistribution of the carrier density, coupling the y and 
z degrees of freedom, reducing the spatial symmetry, and 
inducing orbital hybridization, as discussed in Sec. lIII Al 
From this point of view, increasing is not equivalent to 
decreasing Lf,, since it induces GS configurations different 
from those occurring at large Lb- Moreover, as we show 
below, the high- i^y phases are brought about by genuine 
correlation effects, and in fact they are missed by the 
SSD approximation. 

Let us consider the spin phase at high Bn and low B± 
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in Fig. 13 having the Bu orbital doubly occupied in its 
most-weighted SD. At high the orbitals Ag, are 
almost degenerate (see Fig. EJ. Therefore, we would ex- 
pect the two configurations with two electrons either in 
the Ag or in the i?„ level to occur with like probabili- 
ties. However, intra-orbital Coulomb repulsion is lower 
for the more delocalized B^ state. This is shown for B3 
in Fig. IS^b) , where we report the behavior of Uaa [de- 
fined in Eq. ifT^ ]. where a takes the value Ag and 
respectively. The modulation of UAgAg and Ub^b^ with 
B\\ should not come as a surprise, since the SP orbitals 
are strongly affected by the field (Fig. ^ , with the ap- 
pearance of nodal surfaces even in the lowest-energy SP 
orbital Ag. These integrals allow to estimate the value 
of the Coulomb repulsion between two electrons both sit- 
ting in the Ag or in the i?„ orbital, as a function of the 
in-plane field B\\ . In Fig. |HIa) we show the ground- and 
lowest excited-state energies for B3 vs. B\\. The larger 
values of UAgAg at intermediate fields partly explains the 
occurrence of the SD with two electrons sitting on the 
Bu orbital. On the other hand, such GS is not observed 
in a SSD picture, i.e., neglecting correlation, as shown 
in Fig. [71 (note that the B\\ field range is larger in Fig. 
than in Fig. ITJ. At high values of B\\, UAgAg — Ub^b^ 
and the state with two particles on the B^ level becomes 
almost degenerate with the one with two particles on the 
Ag level [Fig. |S{a)]. The results for Bl and B2 samples 
(not shown) are similar. 

The reduced Coulomb repulsion between electrons on 
the Bu orbital is a key concept to understand the stability 
of the unpolarized phase S = 1/2. Indeed, the SD with 
the Ag orbital singly occupied and the i?„ orbital doubly 
occupied is always present in the GS when applying 
(even though it is not indicated in the graphs unless it 
is the dominating one), together with other important 
configurations. The reason is that correlation, namely 
the mixing of different SDs, allows the system to lower 
its energy. At high B\\ , the SD mentioned above becomes 
the dominating one for all considered samples, although 
close in energy to the SD with the Ag orbital doubly 
occupied. As in other cases, in a tilted magnetic fields 
the typical SDs appearing in the GS expansion, in the 
central region of the diagrams, represent the competition 
between vertical correlation, driven by _B|| , and in-plane 
correlation, driven by B^. 

For = 4 the situation is more complicated due to the 
larger number of possible phases, and the strong depen- 
dence on the inter-dot coupling strength. SP, SSD, and 
FCI predictions are presented separately in Figs. [HI El 
and ^2 respectively. The predictions of the SP scheme 
in Fig. (51 if compared with those of FCI of Fig. ^] are 
clearly inadequate to describe the physics of the system: 
only two non-interacting singlet phases are stable, differ- 
ing for the double occupancy of either the Au or the Bu 
orbital. 

As expected, the FCI calculation (Fig. Ill|l predicts 
the familiar MDD state formation through successive 
occupation of p and d orbitals; this is also predicted. 
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FIG. 9: A'^=4 stability diagram in the {B^\,B±^ plane for 
Bl, B2, and B3, as obtained within the SP non-interacting 
scheme. 




FIG. 10: 7V=4 stability diagram in the (_B||,_Bx) plane for 
B2, calculated within the SSD scheme. 

even though approximately, by the SSD approximation 
(Fig. EU- On the other hand there are many differences 
between SSD (Fig. ^ and FCI (Fig. El) predictions 
when is also increased. One major difference is the 
enhanced stability of the = phase in FCI calculations. 
This phase involves the SD made of doubly occupied Ag 
and Bu orbitals, which is the dominant one in a large field 
range for all samples. Analogously to the A = 3 case, 
the system is allowed to lower its energy through occu- 
pation of the Bu level as B\\ is increased. The true GS, 
however, turns out to be a non trivial linear combination 
of several other SDs, which suggests that correlations are 
dominant here. Indeed, by neglecting correlation effects, 
like in the SSD approach, the only possible energy gain 
is obtained by aligning the spins in a S* = 1 phase, lead- 
ing to the wrong SSD prediction that the S = 1 phase is 
always stable at small B^, as shown in Fig. 1101 

Such partially polarized phases, instead, are typically 
stable only in small ranges of the field, and their stability 
also depend on the sample parameters, as we show be- 
low. For N — A: the S = \ phase at i? « is stable only 
in a small field range, and its stability is reduced if the 
inter-dot coupling is reduced; in sample B3 this phase is 
not stable, no matter how small the field is fFig. lll|l . In- 
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FIG. 11; N — 4: stability diagram in the {B\\,B±) plane for 
Bl, B2, and B3, as predicted by the FCI calculations. Gray 
triangles: S — 2 GS; gray squares: S = 1 GS; black circles: 
S — GS. Solid lines are guides to the eye. Insets: main 
components of the GS wave functions (see text). 
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FIG. 12: N — 5 stability diagram in the {B^^,B±) plane for 
Bl, B2, and B3, as obtained within the SP non-interacting 
scheme. 




deed, this phase requires occupation of two S p orbitals, 
and it results from the payoff between the gain in ex- 
change energy and the kinetic energy cost. Therefore, 
this configuration becomes rapidly unfavorable with re- 
spect to the S = phase having the i3„ orbital doubly 
occupied, as increases. In fact, the stability range 
of the S — 1 phase decreases as the zero field tunnel- 
ing energy Asas decreases with respect to the s — p en- 
ergy splitting hioo. Note that, according to the SSD data 
(Fig. ^J, the 5' = 1 phase is given by only one SD, in- 
stead of the FCI superposition of two main configurations 
with coefficients dependent on the value of the field. This 
fact, which indicates the progressive filling of the level 
while increasing the field, is neglected by the less accu- 
rate SSD picture, and it is described by SSD model in 
terms of a transition between the two SDs mentioned 
above. We note in Fig. ^2 that a S" = 1 phase sets in 
also between the unpolarized {S = 0) and the completely 
polarized (5 = 2) spin case. This phase, characterized 
by a non trivial combination of configurations involving 
both S and AS states, loses its stability as the inter-dot 
barrier is increased. 

Similar considerations apply to the N ~ 5 case (see 
Figs, dill and m for SP, SSD, and FCI schemes, re- 
spectively). Typically, several transitions occur in the 
GS within the same spin configuration. For example, 
along the Bii axis of Fig. E| several configurations oc- 



FIG. 13: N — 5 stability diagram in the (B^^,B±) plane for 
B2, calculated within the SSD scheme. 



cur, reflecting the crossing between p levels in the S and 
AS mini-bands (see Fig. |2J|. As the barrier Lb is de- 
creased, the splitting between s and p orbitals becomes 
of the same order of magnitude of AsaSi and the num- 
ber of transitions in the GS increases, too. Also in this 
case a SSD approach is inappropriate. As in the = 4 
case, we have an intermediate phase {S = 3/2) between 
the unpolarized {S = 1/2) and the completely polarized 
{S — 5/2) ones, which is only stable in the intermediate 
coupling regime; this phase is absent in the SSD results 

(Fig. eg. 

The above results were obtained by keeping huo fixed 
with respect to N. In single-electron tunneling experi- 
ments, however, Tiujq changes with the gate voltage, in 
order to inject (extract) electrons into (from) the AM.— 
A relation linking luq with N may be derived by impos- 
ing that the electron density is approximately constant 
during the charging processes, ujq oc N~^/'^ (Refs. "O^T^. 
This allows us to qualitatively understand what changes 
the AM phases undergo if we try to mimic the depen- 
dence of Wo on N. Indeed, if Uujq decreases as N in- 
creases, the ratioes of this energy to the other relevant 
energy scales — namely AsaSi ?i'^c, and Coulomb matrix 
elements, respectively — decrease as well. This implies 
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FIG. 14: N — 5 stability diagram in the {B^^,B±) plane for 
Bl, B2, and B3, as predicted by the FCI calculations. Gray 
triangles; S — 5/2 GS; gray squares: S = 3/2 GS; black 
circles: S = 1/2 GS. Solid lines are guides to the eye. Insets: 
main components of the GS wave functions (see text). 



that: (i) the AM phase boundaries in the (-By, B±) plane 
move in order to "squeeze" the stability regions of "ionic" 
phases (ii) the transition to the MDD phase along the B±^ 
axis is faster (iii) correlation effects are stronger. 

Finally, we note that the present results for iV = 4, 5 
at = are in agreement with the ones shown by 
Rontani et al. dealing with similar samplesfi^ except for 
little discrepancies due to different parameters used in 
the calculations, in particular the band-offset Vo and the 




lateral confinement energy Tiloq. 



IV. CONCLUSIONS 

We have investigated the relative stability of quantum 
phases of interacting electrons in AMs under a magnetic 
field. Such phases can be identified, e.g., in transport 
experiments in the Coulomb blockade regime, where the 
addition energy is measured.^ A very sensitive test of few- 
electron phases is non-linear transport through charged 
QDs, where a few excited states in a small energy window 
are accessediiiiSiii In addition to predict the stability re- 
gions of the different spin phases, this work confirms that 
correlation effects play a fundamental role in determining 
the physics of the system, and a simple picture neglect- 
ing these effects could fail in interpreting the relevant 
experiments. 

In particular, we have compared FCI predictions with 
(i) a SP picture, neglecting all interactions, (ii) a SSD 
scheme, where only direct and exchange interactions be- 
tween charge carriers are taken into account. These com- 
parisons allowed us to identify the effects of correlations, 
particularly those induced in the coupled structure by 
an in-plane magnetic field, in analogy to what happens 
in single QDs with an applied vertical field. In many 
respects, only the exact diagonalization of the few-body 
Hamiltonian is able to provide a correct description of 
correlated states. 
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